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Abstract 

We present a theoretical description of excitons and positively and negatively charged trions in 
"giant" CdSe/CdS core-shell nanocrystals (NCs). The developed theory provides the parameters 
describing the fine structure of excitons in CdSe/CdS core/thick shell NCs as a function of the 
CdSe/CdS conduction band offset and the CdSe core radius. We have also developed a general 
theory describing the fine structure of positively charged trions created in semiconductor NCs 
with a degenerate valence band. The calculations take into account the complex structure of the 
CdSe valence band and inter-particle Coulomb and exchange interaction. Presented in this paper 
are the CdSe core size and CdSe/CdS conduction band offset dependences (i) of the positively 
charged trion fine structure, (ii) of the binding energy of the negatively charged trion, and (iii) 
of the radiative decay time for excitons and trions. The results of theoretical calculations are in 
qualitative agreement with available experimental data. 



PACS numbers: 73.22.-f, 78.67.Bf, 81.05.Dz 



I. INTRODUCTION 



Growing attention to CdSe/CdS core/thick shell nanocrystals (NCs) initially reported 
by two groups^ is stimulated by their superior optical properties superior to any other NC 
hetero-structures prepared up to now. The photoluminescence is never completely quenched 
in these NCs^ and the blinking is almost completely suppressed^. At low temperatures these 
structures demonstrate suppression of non-radiative Auger recombination^ and almost 100% 
photoluminescence quantum yield^. As a result these NC structures demonstrate a very low 
threshold for optically pumped lasing£. 

The origin of such outstanding optical properties of the CdSe/CdS core/thick shell nano- 
structures (or giant CdSe/CdS NCs^) is not completely understood, but one thing is clear. 
The very thick shell always prevents carrier escape to the NC, no matter which mechanism 
is responsible for this escape: direct quantum mechanical tunneling, thermo-ionization of 
the NCs, or Auger auto-ionization of the NCs, where two electron hole pairs were created 
simultaneously^. The prevention of the NC ionization could explain the blinking suppression; 
however, it does not explain the suppression of non-radiative Auger recombination of charged 
excitons, biexcitonsi^&S and multi-excitons^ 1 ^. 

The alloying of the CdSe/CdS interface in giant CdSe/CdS NCs could be a reason for 
Auger recombination suppression, as was suggested in several papers of the Los Alamos 
group£~-^. This explanation relies on a theoretical prediction that softening of the 
carrier confined potential could reduce the Auger recombination rate by three orders of 
magnitudeii^. Indeed, growth of the thick shell in the giant CdSe/CdS core/shell NCs 
lasted for a week^~-^. The slow growth of the shell at temperatures above room temper- 
ature suggests that the strong interface diffusion takes place during the NC growth. The 
diffusion should lead to alloying adjacent to the surface area and softening of the confine- 
ment potential. This alloying assumption was confirmed recently by the demonstration of 
mixed phonon replicas in the fluorescence line narrowing (FLN) of CdSe/CdS NCs with a 
thick shell^. 

This argument, however, cannot be used to explain the suppression of the Auger re- 
combination observed by Dubertret groupie. The shell of their NCs is grown very fast 
and that practically excludes any inter-diffusion. In contrast, the fast growth of the lattice 
mismatched CdSe/CdS core/shell NCs should obviously results in some intrinsic strains. 
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Indeed, the samples exhibit a phonon spectrum of a pure CdSe NC, the frequencies of which 
are shifted from bulk CdSe frequencies due to the compressive strain within the NO^. 

An understanding of the unusual and potentially useful properties of the giant CdSe/CdS 
core/shell NCs requires a theoretical description of the fine structure of the band edge ex- 
citons, and trions. Up to now, the energy spectra of these NC quantum dots were studied 
only within simplified parabolic band models^ 1 ^ 1 ^, which do not take into account the com- 
plex structure of the valence band and the complicated character of inter-particle Coulomb 
interactions. The theoretical description of the energy spectra of the exciton and exciton 
complexes in these NCs also has some objective complications, because CdSe/CdS hetero- 
structures are not well characterized. We know that holes are always localized in CdSe, and 
that the valence band offset is no less than 500 meV. However, we don't know the sign of 
the band-offset of the conduction band nor its magnitude. A small value of the CdSe/CdS 
conduction band offset suggests its complex behavior, because the band offset sign and 
magnitude should be sensitive to the strain and temperature^. 

These qualitative data on the band offset suggest that an exciton created in CdSe/CdS 
NCs with a thick CdS shell should look like a donor center, where a hole is strongly localized 
in the CdSe core while an electron is attracted to this charged core by the Coulomb potential. 
If CdSe/CdS is the type II hetero-structure, the core potential in the NC pushes the electron 
into the CdS shell. In the opposite case, if CdSe/CdS is the type I heterostucture,the core 
creates an additional well, which increases the electron localization around the hole. The 
temperature dependence of the CdSe and CdS energy gaps^ suggests that a decrease of 
temperature should favor type I. As a result, a decrease of the temperature increases the 
exciton binding energy and, consequently, shortens the radiative decay time of excitons, 
which is controlled by the electron-hole wave function overlap, and increases the dark-bright 
exciton splitting.— 

Due to the very strong confinement of holes, the positively charged trion which was 
observed in the CdSe/CdS NCs^ appears as a doubly-charged donor. The larger Coulomb 
attraction of the electron to the core should significantly increase the electron binding energy 
and shorten the radiative decay time. The negatively charged exciton should look like a 
negatively charged donor center D ~~ . The binding energy of the second electron in the 
D ~ center is usually quite small on the order of ~ 0.055 of the donor binding energy. 
This energy will be higher, however, in redthe type I CdSe/CdS core/shell structure, due 
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to an additional attractive potential (created by the CdSe core) for an electron. Finally, a 
biexciton in these structures looks like the helium atom, with two holes localized in the core 
and two electrons with opposite spins occupying the same orbit. 

In this paper, we have calculated the binding energy and the fine structure of the exciton 
in CdSe/CdS NCs with a very thick shell as a function of the CdSe/CdS conduction band 
offset and the core radius. The calculations takes into account the complex structure of 
the CdSe valence band and inter-particle Coulomb and exchange interaction. We have 
also developed a general theory describing the fine structures and radiative decay time of 
positively and negatively charged trions. Using this theory we have calculated optical spectra 
of trions in giant CdSe/CdS core/shell NCs as a function of conduction-band offset and core 
radius. 

This paper is organized as follows: In Section II and III we describe the energy spectra 
and wave functions of one and two holes localized in the CdSe core. In Section IV we 
discuss the energy and the wave function of an electron in a central symmetrical potential 
created by localized holes and the conduction band offset. In Section V we derive the wave 
functions and the energy of excitons and positively charged trions. In Section VI we discuss 
the fine structure of these excitons and trions. In Section VII we calculate the binding 
energy of the negatively-charged trions. In Section VIII we consider the radiative decay of 
all exciton complexes. Finally, we discuss obtained results and compare them with available 
experimental data in Section IX. 

II. HOLE ENERGY SPECTRA AND WAVE FUNCTIONS 

In CdSe/CdS core shell nanocrystals, the valence band offset is substantially large and the 
holes are strongly confined in the CdSe core. For the heavy ground state holes this potential 
can be considered as infinitely high in the first approximation. The first quantum-size level 
of holes in a spherical NC of a semiconductor with the degenerate r 8 valence subband 
is a IS3/2 state.— This state has total angular momentum j = 3/2 and it is fourfold 
degenerate with respect to its projection M = 3/2, 1/2, —1/2, —3/2 on the z axis.—. The 
wave functions of these four states can be written as^ 

^ = 2EH)^« E ( l W«m> (!) 

1=0,2 m+ti=M \ m V ~ M J 
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where k n p z ) are Wigner 3j-symbols, Y^ m (9, (/)) are spherical harmonics defined in Refil9 
and u v (// = ±1/2, ±3/2) are the Bloch functions of the fourfold degenerate valence band 
r 8 — . The radial functions Rq{t) and R2{f) corresponding to an impenetrable barrier at 



the CdSe/CdS interface with radius a can be written in explicit form following Refs. 
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where jo an d H are spherical Bessel functions, the constant A is determined by the nor- 
malization condition J Q a drr 2 [R^{r) + R 2 {r)] = 1 and j3 = mi/m^ is the ratio of light to 
heavy hole effective masses: mi and , respectively. The dimensionless parameter (p is 
a first root of the equation jo(v 9 )i2(\/^^) + j2{<p)jo(V~P l P) = 5 anc l describes the energy of 
the lowest hole level^: £as 3/2 = h 2 ^ 2 / (2m h a 2 ) . 

The assymetry effects, which include the intrinsic asymmetry of the hexagonal lattice 
structure of the crystal's field,— and the nonspherical shape of the NC^ 3 ., split the four-fold 



degenerate hole state into two two-fold degenerate states with \j z 
respectively: 



AOS, a) 



E lSm , Su = E lS3/2 - -^{j 2 z - 5/4) 



3/2 and \j z \ = 1/2 



(3) 



where A(/3, a) = A int (/3)+A sh (/3, a) where A int and A sh are the intrinsic and size dependent 



shape contributions which can be found, for example, in Ref. 
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III. TWO HOLE ENERGY SPECTRA AND WAVE FUNCTIONS 

Before considering the positively charged excitons, which can be observed in NCs with 
an extra hole, let us consider the energy spectra of two holes localized in the CdSe core. 
The application of the Pauli exclusion principle to the two holes occupying the 15*3/2 level 
is nontrivial. Generally, the two holes with momentum \ji\ = \j 2 \ = 3/2 according to 
the momentum summation rule J — ji + 32 could form four states with the total angular 
momentum J = 3, 2, 1 and . The wave functions of these states can generally be written 
as 

(3/23/2 J \ 
U^(n)^ 2 (r 2 ), (4) 
M l M 2 -J Z J 
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where J z is the projection of the total momentum J on the z axis. The Pauli hole 
permutation requirement applied to the wave function described by Eq. (j3J), however, allows 
only the two nontrivial solutions. As a result, the two holes occupying the 15*3/2 level can 
only be in a 5-fold degenerate state with total momentum J = 2 and a state with J = Mi2k 
Hole-hole exchange interaction splits the ground biexciton into the two states with total 
momentum J = 2 and J = . The straightforward calculation of the h — h Coulomb 
interaction with the functions from Eq. (jl]) gives the energy of the corresponding levels: 

Ejtf) = E 2h {(5) + A cxch (/3) (\ - J(J 6 +1) ) (5) 

where E 2 h((3) = 2Exs 3/2 + E^ is the band-edge energy of the two- hole state taking into 

account the direct Coulomb repulsion only, A cxch (/3) = E (f3) — E 2 (f3) is the exchange 

splitting and can be written in the following form 2 ^ 2 ^ 

e 2 32 f f r 2 
A exc h(/3) = ~-n rlr 2 2 dr x dr 2 -f R^R^R^R^) , (6) 
k 25 J J r% 

where k is the dielectric constant of the semiconductor and r< = min{ri,r 2 } , r> = 

max{ri,r 2 } . The direct Coulomb interaction energy is given by 

E 2 ^m = ^- J j rlrldr 1 dr 2 ^{Rl{r 1 ) + Ri{rMRl{r2)+ (7) 

In the spherical NCs, the ground two hole state has total momentum J = 2 .— In the case of 
the CdSe/CdS structure with impenetrable interface, when the wave functions are described 
by Eq. 02]), the exchange splitting and Coulomb energy can be written as 

A CX ch(/3) = —7(0) , = ~XW) ■ (8) 

KCL KCL 

In CdSe NCs, where (3 w 0.28, 7(0.28) w 0.033^4 and x(0.28) w 1.92. 

The NC asymmetry, which lifts the degeneracy of the IS3/2 hole state in Eq.([3]), causes 
further splitting of the two hole states. The corresponding perturbation can be written as 2 ^ 

Hi h = HLHn) + H^(r 2 ) = - 5/4) - ^^Qi - 5/4) , (9) 

where j± z and j 2z are the operators of hole momentum projection acting on the first and 



second hole, respectively. The resulting two-hole states are described in Ref. [2J. The ground 
two-hole state has angular momentum projection of J z = and energy 



A exc h(/3, a) I ( A exch (/3, a) x 



The first excited state has J z = ±1, and ± 2 and energy Ei , which is not affected by 
perturbations connected with NC asymmetry. Finally, the upper state has J z = and the 
energy 

77+ 77 fa \ Acxch(/3,a) // A cxch (/3,a)\ 2 

E£ = E (P,a) + WI - I + A(/3,a) 2 , (11) 

Two-hole wave functions $Q h,± (for the energies ) and $^ ±2 (for the fourfold degen- 
erate intermediate level _E 2 ) are given by^ 

^(rifl, r M ) = i + S ±) (^M^M - *^ /2 (r M )^ /2 (r, 2 )) + 
(S± - B±) (*? /2 (r M )tt* 1/2 (r M ) - ^i/aCrfci)*^^))] , 
*£(r«,r M ) = ±-^= (^/aCrfcO^Crw) - ^ 1/2 (^1)^3/2(^2)) , 

$S(r w , r M ) = ±-^= (^±3/2(^)^1/2(^2) - ^1/2(^1)^3/2(^.2)) , (12) 

where 

B ° ±==F vAr B2± = ^ - (Bf )2 ■ (13) 

with D± = A cxch =f v^+W. 

IV. ELECTRON ENERGY SPECTRA AND WAVE FUNCTIONS 

In Section II above we have considered holes in the CdSe/CdS core/shell NCs assuming 
that they are very strongly confined in the CdSe core. In that case their wave functions are 
not modified by the Coulomb interaction with the electron or by the hole-hole interaction, 
and are solely determined by the radius of the CdSe core. On the contrary, the conduction 
band offset in such structures is known to be small and even the sign of the band offset is 
not confirmed. As a result, the band offset itself does not localize (or only weakly localizes) 
an electron in the CdSe core. The attractive Coulomb potential created by the holes plays 
an important role in electron localization at the center of the core/shell NC structures. 

We will show in the next section that the attractive Coulomb potential created by holes 
localized in a spherically symmetrical CdSe core does not have spherical symmetry. In the 
first approximation, however, the ground electron state is described by the wave function of 
S symmetry, which can be written as 

ni/aW = Re (r)F o(e)< 1/2 , (14) 
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where u c s = \S,S Z > are the Bloch wave functions of the conduction band and S z = ±1/2 
is the electron spin projection. The radial wave function R z (r) is described by the radial 
Hamiltonian 

^=-^T)hU^) +v '^ (15) 

m e (r) = mcdSe, V z {r) = V z oul (r) + U os , r <a, 
m e (r) = m C ds, V z {r) = V z oul (r) . r > a 

Here Uqs is the conduction band offset (CBO) ( Uqs > in type II structures and Uos < 
in type I structures), mcdSe and mcas are the electron effective masses in CdSe and CdS, 
correspondingly, and V z oul (r) is the effective Coulomb potential of one (Z = 1) or two 
( Z = 2 ) holes acting on an electron. We use Z = to describe the "free" electron problem 
in the case of the neutral core. 

The radial wave function R z [f) is completely defined by Eq. (fl5|) and the standard 
boundary conditions at the CdSe/CdS interface r = a: 

R z (r)\ r<a = R z (r)\ r>a , 
1 dR z (r) , _ 1 dR z (r) , 

j \r<a — , \r>a ■ K 10 ) 

m CdSe dr m CdS dr 

Later we will need the radial function of the "free" resident electron, which remains in 

NCs after negative trion recombination. Its wave function R° e (r) can be written as 

hi (ina) . . 

-jo {kr) r < a 



Rl(r) = A{ jo(ka) J ~ (17) 



h\ (inr) r > a 



where A is the normalization constant, Hq is the spherical Hankel function,- 
yj— 2m CdSc E® j h , k = ^/2m CdS (E® + Uos)/h and E Q e is the electron ground state en- 
ergy level calculated from the bottom of the CdS conduction band. This energy is described 
by the smallest solution of the following equation: 

1 — ka cot (ka) — \i (1 + na) (18) 

where /i = rft CdSe /m C ds • If is important to note that even if CdSe/CdS core/thick shell NCs 
are Type I structures, the resident electron is not localized in the CdSe core with small radius 
a and is spread over the entire CdS shell. The electron localization occurs only if the power of 
a quantum well w = ay/2mcdScUos/h is larger than , which is a solution of the following 
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equation: cot(wQ r ) = (1 — h)/wq . For /x = 1 the electron derealization occurs when 
= 7r/2 , the condition which is identical to the standard expression for the critical depth 
of the 3D potential allowing the electron localization, Uos = ir 2 h 2 / (8mcdSe a2 ) ■— However, 
the large electron effective mass in CdS resulting in \i = 0.634 allows the localization in the 
CdSe core with w smaller than 7r/2: it/2 > w > w^ 1 = 1.296 . 

V. DONOR-LIKE EXCITON AND POSITIVELY CHARGED TRION 

The electron wave function is strongly modified by a long range Coulomb potential created 
by a hole or holes in the donor- like exciton or the positively charged trion in CdSe/CdS NCs 
with a small CBO. Electrons are bound to the CdSe core by a strong Coulomb potential of 
a single hole or two holes, Z = 1, 2 correspondingly, which are strongly confined by the core 
potential, similar to an electron in singly or doubly charged donors. This potential always 
leads to some binding of electrons to the CdSe core, even in type II CdSe/CdS core/shell 
structures. The attraction potential of holes is described by the distribution of the hole 
positive charge within the core of the radius a . 

A. Hole charge distribution and Coulomb potential 

The straightforward use of the hole wave function definition in Eq. ([I]) gives the charge 
distribution of the single hole state with M = ±3/2 and M = ±1/2 correspondingly: 




Using the definition of spherical harmonics, Yi >m , frorn^ we arrive to 

p M (r) = -jL [R 2 (r) ±2(|M| - l)R (r)R 2 (r)(l - 3 cos 2 #) + R 2 2 (r)} . (20) 

The charge distribution of two isotropic holes with A = , localized in CdSe core, p 2 j l j z {v) 
can be obtained from the two-hole wave functions $j h j 2 defined in Eq. (j3j): 

f%j.{r) = J dh hl dh h2 \^f Jz (r hl ,r h2 )\ 2 [S(r-r hl ) + S(r - r h2 )} (21) 
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To obtain the hole charge distribution for A ^ one has to use the two hole wave functions 
described in Eq. (|12p . The calculations show that for the four-time degenerate hole state 
E 2 , which is characterized by the angular momentum projections J z = ±1, ±2 , the charge 
distribution is spherical and is described as 

P±U(r) = ^ + Rl(r)} (22) 

For the upper, " Eq " , and lower, " Eq " , two hole states with J z = similar calculations 
result 

ftf*(r) = A [R 2 (r) + Rl{r) + 2B±B±R {r)R 2 {r){l - 3 cos 2 9)] (23) 
In the limit A 3> A exc /i Eq. (J23]) is simplified significantly: 

pf^(r) = |- [i? 2 (r) +i? 2 (r) Tj R (r)i? 2 (r)(l -3 cos 2 ^)] . (24) 

Comparison of Eq. (|24|) with Eq. (1201 shows clearly that in the case A 3> A exc h , the lowest 
level is created from two heavy holes, while the uppermost level is created by two light holes. 
Similar comparison of Eq. (1221) with Eq. (|20|) shows that the middle level is always created 
by one light and one heavy hole. 

One can show that only the spherical part of the hole charge distribution contributes to 
the ground S electron state in first order perturbation theory. For the single ( Z = 1 ) and 
double ( Z — 2 ) hole states this distribution can be written as 

Po(r) = ^[Rl(r) + Rl(r)] (25) 
The Coulomb potential created by the hole spherical charge distribution can be written 

as 

k Jq \r — r'\ 

where k is the dielectric constant, which is considered to be the same in both CdSe and 

CdS. Equation fl26|) can be rewritten as: 

Ze 2 4vr [1 
k r 

This potential outside the CdSe core, r > a , always has a standard Coulomb orm, V(r) = 
Ze 2 /(m") , due to the wave function normalization condition, f Q a drr 2 [R^{r) + R 2 {r)] = 1 . 
Inside the core V'(r)\ r=0 = 0, and the potential depth V(0) = —{Ze^ATi/n) drrp^r) 
depends only on the ratio of the light- to heavy-hole effective masses. In CdSe where 
(3 = 0.28 ; V(0) = 2.684Ze 2 / (no) . 
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v&Jr) 



dr'r' 2 p z (r')+ I dr'r'p^r') 



n 



(27) 



B. Electron wave function and binding energy 

We will consider only CdSe/CdS core/shell NCs with a very thick shell because these 
NCs show the most interesting and useful properties. From now on all calculations will be 
conducted for heterostructure with unlimited CdS shell thickness; later we will discuss the 
limitations of our approach. 

The ground state of a weakly bound electron in the donor-like exciton or positively 
charged trion created in CdSe/CdS is described by the Hamiltonian of Eq. (II 5p with 
Coulomb potential, Vq ou1 , described by Eq. (1271) . Introducing dimensionless distance from 
the NC center, x = r/a, and dimensionless energies E = E j \h 2 /2mcdSe a2 ) % we obtain the 
dimensionless Hamiltonian, 



use V(x) = 2Vc ou \{r) / (Ze 2 / no) , Uos = Uos/(h 2 /2mcdSeQ> 2 ) ■ All energies are calculated 
from the bottom of the CdS conduction band, and the confinement potential Uqq created 
by the conduction band offset is negative in type I, and positive in type II, CdSe/CdS 
hetero-structures. We also introduce the unitless parameter rjo = clb/cl that characterizes 
the strength of Coulomb potential created by holes, where a# = h 2 n/ (mcdSe^ 2 ) is the Bohr 
radius. 

A central parabolic part of the hole potential described by Eq. (1271) becomes important 
if Z/rj is large. The wave function of the IS ground state should be close to the wave 
function of the harmonic oscillator. In the opposite case of small Z/t]q , the core potential 
is weak and the electron moves mainly in the Coulomb-like potential at large distances from 
the core. To describe both above-mentioned limits, we select a trial radial function in the 
following form: 



where a and ( are dimensionless variation parameters, which will be found using the 
variational principle, and constants A ,B , and D , which are completely determined by the 




(28) 



for the radial component of the electron wave function, i?f 



(x) . In Hamiltonian H r (x) we 




(29) 
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interface boundary and normalization conditions. The standard boundary conditions at the 
CdSe/CdS interface (x — 1) described by Eq. ffTBT) gives 



Rt{a,C;x) 




1+x2 (2C-«V) 



Z( n f.^- ' I 2 (« 2 - 1) C + «VJ 



2a 2 C 

x > 1 



x < 1 

(30) 



2 (a 2 - l)C + « 2 /i 

where B can be found from the normalization condition J °° dxx 2 R 2 e (a ) £; x) = 1 (see 
Appendix). 

Using Eq. fl28|) we can find the expectation value for the energy, E z , of one electron: 

oo 

(E z ) e = J R z (x)H z (x)R z (x)x 2 dx (31) 
o 

To simplify the calculation we approximate the potential, V(x) which is inside the core 

( x < 1 ) by the polynomial V(x) ~ p + p 2 x 2 + p 3 x 3 + p 4 x 4 + p 5 x 5 , where p = —5.368 

and pi — were selected to satisfy the exact potential properties at x = . Outside the 

core (x > 1) the potential is given by V = —2/x according to Eq. (l27|) . The potential 

and the potential derivatives should be continuous at x = 1 , which allows us to determine 

Pi = 14.84257478-3p 2 -2p 3 and p 5 = -11. 4740598 + 2p 2 +p 3 . Next, we vary the coefficients 

P2 and p-s to fit the polynomial to the function V (x) , the form of which depends only on the 

parameter /3 . The best fit of V (x) for = 0.28 was achieved using p 2 = 10.6559838 and 

P3 = —7.111613 . Figure [TJ shows the dimensionless potential V (x) . The COB is introduced 

via [Jos One can see that when the band offset is positive, the electron is pushed away 

from the CdSe core. The negative band offset added to the Coulomb attractive potential 

stimulates electron localization. The second hole localized in the CdSe core just increases 

the magnitude of the attractive potential by a factor of two. 

Using this approximation we were able to develop an analytical expression for the electron 

binding energy (E z ) e (a,() , via parameters a and ( and find the absolute minimum of 

the total energy and a corresponding pair of ( and a for given /i , r] , Z , and the band 

offset £7off (see online Supplemental Material). 

Figure |2] shows the binding energy of an electron for various core radii and band offsets, 

calculated for the cases when the CdSe core is charged with one ( Z = 1 ) and two ( Z = 2 ) 

holes. The energy is calculated from the bottom of the of the CdS conduction band. The 

size dependences are calculted for a CdSe core with a 2.5 nm radius and shown for five band 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Dimensionless radial distance, r/a 

FIG. 1. Dimensionless Coulomb potential V (x) created by one hole {Z = 1). Insert shows 
schematically the total potential created by the sum of the Coulomb potential acting on the electron 
and the conduction band offset. 

offsets: u = 50, 0, —50, —100, —200 meV. One can see that in NCs with small negative or 
positive band-offset, this energy decreases with a radius increase, while in the case of -100 
and -200 meV, it increases following the bottom of the conduction band. The decrease of 
the binding energy with radius increase in NCs with small negative or positive band-offset is 
connected with the fact that electron motion in such NCs is controlled by Coulomb potential, 
which decreases with NC size. 

Figure [2] shows the dependence of the electron binding energy on the band offset, Uos , 
in NCs with core radius 1.5, 2.5, and 3.5 nm radius for the single and doubly charged CdSe 
core. The band offset depends on temperature^, and the shift of the exciton line can be 
observed in temperature dependence. 
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FIG. 2. Binding energy of an electron in the Coulomb potential created by one and two holes 
localized in the CdSe core calculated for various core radii and conduction band offsets. The 
dependence of the electron binding energy on the CdSe core radius, a , for one hole (Panel a) and 
two holes (Panel c). The dependence of the electron binding energy on the conduction band offset 
Uos for one (Panel b) and two holes (Panel d). 

VI. THE FINE STRUCTURE OF EXCITONS AND POSITIVELY CHARGED 
TRIONS 

The fine structure of the excitons and trions is strongly affected by the electron-hole 
exchange interaction, which can be written in spherical approximation as 

2 

# CX ch = ~-al5(r e - r h )e cxch (er ■ j) , (32) 

where, ao is the lattice constant, 5 cxc h is the strength of the exchange interaction (in energy 
units), er are the Pauli matrices, and j = j x ,j y ,jz are three matrices of the momentum 
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3/2. The strength of the exchange interaction in cubic and hexagonal semiconductors can 
be determined from the bulk exciton fine structure. In CdSe this constant extracted from 
the bulk exciton splitting e cxch = 450 meV.— 

Averaging Eq. f[3"2"j) over electron and hole wave functions in the basis of the 8-fold degen- 
erate exciton stated, one can arrive at the following expression for the spin-spin exchange 
Hamiltonian between an electron and a hole: 

HU = -Vi(° ■ 3) i (33) 

where the exchange parameter rjz depends strongly on the number of holes in CdSe core: 

Vz = ^e exch jyrr 2 [R^(r)} 2 [R 2 (r) + 0.2Bi(r)], (34) 

due to the strong effect of the core charge on the electron radial function, -Rf (r) . The radial 
functions i?f (r) and -Ro,2( r ) are defined in Eqs. fl30|) and ([2]), respectively. Integration in 
Eq. ( 134]) goes only to the core radius a because in our approximation the hole wave function 
Ro^fr) vanishes at the CdSe/CdS interface. One can see from Eq. ( 1341) . that if an electron 
is only weakly localized in the CdSe core, then the parameter rj is small. This parameter is 
proportional to (a /aQ 3 <C 1 where a( is a localization radius of the electron. 

In Fig. Qa we show the dependence of the exchange energy t]i of the exciton as a 
function of the conduction band offset for several CdSe core radii. One can see that increase 
of £/off leads to an increase of electron localization in the core always increases rji . Figure 
([3])b we show the dependence of exciton exchange energy rji as a function of the CdSe core 
radius for several NCs band offsets. For all reasonable set of radii, rji increases with the 
decrease of the NC radius. 

The ground exciton state in spherical NC of zinc-blende semiconductors is characterized 

by the total exciton momentum F = j + S . For the ground 5-fold degenerate Dark exciton 
state F = 2 , and for the excited 3-fold degenerate Bright exciton state F = 1 . In cubic 
semiconductor NCs with nonspherical shape and NCs with wurzite lattice structure these 
states are split into 5 excitons, each of them characterized by the projection of the total 
mometum on the hexagonal axis F z .— The respective wave functions are constructed from 
the hole functions of Eq. ([I]) and electron functions of Eq. (jl4j) . with the radial functions 
i?f =1 for the electron moving in the Coulomb potential created by one hole localized 

in the CdSe core. 
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FIG. 3. Dependence of the exciton exchange parameter r/i on the conduction band offset, t/off (a), 
and CdSe core radius, a (b). Insets show the same dependences for the trion exchange parameter 
V2 ■ 

The band edge exciton fine structure is completely characterized by two parameters rjx 
and A(a, /3) defined above. In spherical cubic NCs, the splitting between Bright and Dark 
excitons is given by E 2 — Ei = Arji . In the case corresponding to A(a, 0) 3> 771 , the splitting 
between the ground Dark exciton with F z = ±2 and the first excited Bright exciton exciton 
with F z — ±1 is 3r/i . The last case is most probably realized in presently grown CdSe/CdS 
NCs.— In the more general case, corresponding to A(a,B) ~ rji , the exciton fine structure 
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is described by expressions which can be found in Refs. |22|, [24|, and 

The operator for the electron-hole exchange interaction in the positively charged trion 
can be written as 

2 

#cxch(Ti) + H cxch (r 2 ) = --a^cxch [6(r e - r x ){cr ■ ji) + 5(r e - r 2 )(cr ■ j 2 )} (35) 

where r\ , r 2 and ji , j 2 are coordinates and spin matrices of two holes, respectively. The 
averaging of Eq. (135]) over radial components of the two hole and electron wave functions 
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results in the effective electron-hole spin exchange Hamiltonian 

#cxch = -2r] 2 (cr ■ J) , J = j\ + h , (36) 

where 772 is defined in Eq. The dependences of the trion exchange parameter 772 on 

the core radius and the conduction band offset are shown in insets of Fig. fl3]). One can see 
that due to stronger attraction of the electron to the CdSe cores coursed by the second hole, 
the trion's 772 are larger than those for the exciton. 

In the absence of the hole level splitting (A(/3,a) =0), the states of the positively 
charged trion are characterized by the total momentum G — j\ + j 2 + l/2cr = J + S and 
the trion wave functions can be generally found as 

$S, 9 (r e ,r M ,r, 2 ) = (-l) s - 1/2 v / 2GTT £ / J 1/2 J \ ^ 

Jz+S z =g \ Jz S z -g J 

(37) 

where g is the projection of the total momentum G on the z axis, S z = ±1/2 is the 
projection of the electron spin on the z axis, $jj z is the two-hole wave function described 
by Eq. (j4j), and is the electron wave function described by Eq. (fT^j) with the radial 
function R 2 e for the electron moving in the Coulomb potential created by two holes. The 
ground trion state has total momentum G = 5/2 and the energy E5/2 = E 2 — 2772; the 
next one has G = 3/2 with the energy E 3 / 2 = E 2 + 3r] 2 ; and, finally, the upper state 
has G = 1/2 with the energy E\j 2 = Eq . The last trion state is not affected by electron 
hole exchange interaction because the electron interaction with two different holes exactly 
compensate each other. The splitting of the two lowest hole states E 2 coursed by the 
exchange interaction is E 5 / 2 — E 3 / 2 = 5r] 2 . The order of two excited trion states depends 
on the relative strength of electron-hole or hole-hole exchange interactions. All these trion 
states are twofold degenerate with respect of the momentum projection g , and their wave 
functions, g , can be written 

= $0j* ± l/2 , $3/2, ± l/2 = $iI±l^l/2 =F y| $8,^1/2 , 

*3A ± 3/2 = ±)fl (2<^ ±2 *; V2 - S&ini/J , $ 5/2, ± l/2 = + ^^±1/2 , 

*5/2, ± 3/2 = {^±2^1/2 + 2^ ± ^il /2 ) , $^,±5/2 = ^±2^1/2 , (38) 
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When the ground state of holes is split by a crystal field or when the NC shape is 
asymmetric, the related perturbation can be written in the basis of the trion functions of 
Eq. (j38p . The perturbation does not mix the the trion state with different signs of the spin 
projections. For the positive spin projections (g = 5/2,3/2, and 1/2) the perturbation has 
the following form: 





|l/2,l/2> |3/2,l/2) |5/2,l/2) |3/2,3/2) |5/2,3/2) |5/2, 5/2) ^ 
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-\/f A 
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15/2,3/2) 











E 2 - 2r] 2 





|5/2,5/2) 














E 2 - 2r] 2 j 



(39) 



V 

A very similar matrix describes the perturbation created by nonzero A on the states with 
negative angular momentum projections: g = —5/2,-3/2, and -1/2. That matrix can 
be obtained by replacing states \G, g') in Eq. fl39|) by states with \G,—g'), where g' = 
5/2,3/2, and 1/2. The only difference is the sign of the matrix element taken between 



states 1 1/2, — 1/ 2) and 1 3/2, —1/2) , which in the last case is —a/2/5 A . 

One can see from Eq. (1391) . that NC asymmetry does not affect the three trion states 
|5/2, ±5/2) , |5/2, ±3/2) and |3/2, ±3/2) . The energies of these states are 



E5/2 = E 2 - 2r] 2 , E 3/2 = E 2 + 3772 



(40) 



and their wave functions, $3/2 ±3/2 f° r -^3/2 , and $5/2 ±3/2 an< ^ ^5/2 ±5/2 ^ or -^ 5 / 2 are the 
same as in the spherical case. 

The other three trion states, E^ 2 and Ei/ 2 can be found as solutions of the following 
third order equation: 



(E 2 -E- 2r] 2 )(E 2 -E + 3t/ 2 )(£ - E) - A 2 (E 2 - E + m ) = 



(41) 



In the case of small electron-hole and hole-hole exchange interactions < A, A e;EC A 
the approximate solutions of Eq. (1411) yield 



E l/2 



E 2 + A exc /2 ± A , E 1/2 ^E 2 + V2 . 



(42) 
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The wave functions of these trion states are ^!yy2 — ' ^±1/2 ano - ^±1/2 — < ^ ) 2±i^ / ^=i/2 ' 
respectively. 

In the limit of small exchange interaction, the order of the trion levels is independent 
of the relation between A exc and ?] 2 . The fine structure of the positively charged trion 
is, ordered upward in energy, the ground twofold degenerate Ey 2 level, followed by the 
fourfold degenerate E 5 / 2 level, followed by the twofold degenerate Ey2 and E 3 / 2 levels, 
and finally by the upper twofold degenerate Ey 2 level. The lowest and the uppermost trion 
states in the limit of large A are formed from the lowest and uppermost two-hole states 
with a total zero momentum projection of two holes J z = , while all intermediate states 
are formed from the medium fourfold degenerate two-hole level with J z = ±1,±2. The 
medium level of the trion would be eightfold degenerate if one neglected the splitting caused 
by the electron-hole exchange interaction proportional to 772 . 

The limit of small exchange interaction should describe quite well the CdSe/CdS NCs 
which have a core radius ~ 20 A. In these NCs, the hole-hole exchange interaction, A exc ~ 
5 meV and A(/3,a) ~ 20 meV. The calculations of the electron-hole exchange interaction 
in Fig. [3] show that 77 2 is also much smaller than A(j3, a) . 



VII. NEGATIVELY CHARGED TRIONS 



Let us consider now a negatively charged trion formed by the optically excited electron- 
hole pair and resident electron. The total energy of the negatively charged trion can be 
generally described as a sum of the two total energies of electrons interacting with the hole 
described in our case by Eq. ( 13TT) and the energy of the electron repulsion. In the case of 
the small conduction band offsets, the binding energy of the trion is mainly controlled by 
electron-electron and electron-hole Coulomb interactions as is true for the negatively charge 
hydrogen atom. The resulting binding energy in such a case is significantly smaller than the 
binding energy of the exciton because the Coulomb attraction between a second electron 
and the hole is almost exactly compensated by the Coulomb repulsion between electrons. 



The studies of negatively charged hydrogen atoms (see for example Ref. 1291 ) demon- 



strate that the binding state of doubly charged atom exist only for a singlet state of two 
electrons occupying the two radial orbits of S symmetry. This requirement decreases the 
Coulomb repulsion energy between two electrons. That is why we select the two-electron 
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wave functions of the negative trion ^> 2e {ri, 7*2) as the product of the singlet spin wave 
function (antysimmetrical) and symmetrical radial function of two S -symmetry electrons 
in the following form: 

*2e(ri,r 2 ) = -^« /2 (elK_ 1/2 (e2) - n c _ 1/2 (el)< /2 (e2))Foo(e el )r o(e e2 ) J R 2e (r 1 , r 2 ) ,(43) 
where R 2e is defined as 

R2e(r!,r 2 ) = A 2e [Rl(a 1 X\]r 1 / a)Rl{a 2 ,C,2\r 2 / a) + R 1 e (a 1 ,Ci;r 2 /a)R 1 e (a 2 ,C2;ri/a)] (44) 



with Rl(a,(', r AO defined in Eq. (I3"U1) and A 2e defined by the normalization condition 
jo°° jo°° l^2e(ri,r 2 )| 2 r^r 1 r|dr 2 = 1 . 

Using this approach we can write down the total dimensionless energy of two electrons 
in the trion 

/oo r l3 r f jZ r 

i? 2e (r 1 ,r 2 )^ r (r 1 ) J R 2e (r 1 ,r 2 ) 1 -i 1 -^ + ^ oui; (45) 

where the repulsion energy between two electrons E^ oul can be found as 

2e _ e 2 f d 3 ri d s r 2 \R 2e ( ri ,r 2 )\ 2 _ e 2 f , , Rl(r u r 2 ) 
Eco ^ ~TJ 17^T |n-r 2 | " 7 7 r ^ dr ^2 - (46) 

The two-electron trial wave function defined in Eq. (144 j) is described by four fitting param- 
eters, and it generally describes the electron motion on two radial orbits characterized by 
two different sets ai , £1 and a 2 , ( 2 . It is important to note that if a.\ 7^ a 2 and Ci 7^ C2 ; 
the energy of two electrons in Eq. (j45|) is not reduced to the kinetic energy and repulsion 
of uncorrelated electrons occupying two different orbits. The total energy of electrons has a 
contribution from the interference of both orbits and depends on all four parameters. 

One can find the analytical expression of the expectation energy (E) function of 

four parameters in the online SM. This expression developed within Mathematica were used 

to find the binding energy of two electrons and a hole. The trion binding energy, which is 

the difference between the total energy of the trion and the binding energy of the exciton 

(see Fig. |2]), is shown in Fig. H]as a function of of the NC core radius and band offset. 
One can see that for large negative values of Uos , the CdSe/CdS NCs are well-defined 

type-I structures. The strong confinement allows the two electrons to occupy a ground state 

level with two opposite spins leading to a large trion binding energy (see Fig. [5]). Decrease 

of I [/off I decreases the confinement, and the trion binding energy is mainly controlled by 
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FIG. 4. Binding energy of the negatively charged trion in CdSe/CdS core thick/shell NCs. (a) - 
dependence of the binding on the core radius, a , calculated for three band offset [/off =-50, -100, 
and -200 meV. (b) - dependence of the binding on the band offset Uos calculated for three radii 
a = 1.5 , 2.5, and 3.5 nm. Panels c and d show the dependence of the orbit radii £ia and foa, as 
functions of the core radius, a , and of the band offset [/os , respectively. 

the electron-electron and the electron-hole Coulomb interactions. The Coulomb repulsion 
between two electrons forces them to occupy different orbits, which reduces the repulsion 
and leads to a small, but positive binding energy for the trion. The orbits are characterized 
by the two radial wave functions, Ri(a±, £i, r/a) and Rl(at2, (2, r/a) , the spatial extension 
of which is characterized by parameters da and £2° . The radius of these two orbits is 
shown in Fig. |4fc andHJi as a function of core radius and band offset. One can see that one 
of the orbits is mainly localized in the CdSe core while the other explores the CdS shell. 
The radius of the large orbit obviously increases with the decrease of \Uqs\ for any core 
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radius, as one can see in Fig. [4b. The CBO dependence shown in Fig. [4H is more complex, 
however. The radius of the external orbit is smaller in the NC with a large core radius at 
large CBOs, but it becomes larger at small CBOs. This unusual behavior is connected with 
the increased role that the hole Coulomb attraction potential plays in electron localization 
in core/shell NCs with a small CBO. 



VIII. OPTICAL PROPERTIES OF EXCITONS AND TRIONS 
A. Optical transition probabilities of the band edge excitons 

The optical selection rule and the relative transition probabilities of the excitons in 



CdSe/CdS core shell NCs are described by Eqs. (25)-(28) of Ref. |22|, where the square 
of the overlap integral, K , between the electron and hole wave functions should be written 
in more general form 



K, 



ex 



drr 2 R L e (r)R (r) 



2 



(47) 



where R\ is the electron radial function described in Eq. ([29]) . with parameters a and Q 

that maximized the exciton binding energy. In the CdSe/CdS core shell NCs, where the 

electron is not localized in the CdSe core, K ex is a strong function of the band offset of the 

conduction band and of the core radius. In Figs. Eh and [5b we show the dependence of the 

exciton K ex on the core radius and the band offset, respectively. 

The experimentally observed increase of the the radiative decay time by a factor of two 

in the temperature interval from to 150 K— could be connected with a significant 

decrease of the conduction band offset between CdSe and CdS shell layers. The decrease of 

the absolute value of the negative band offset decreases K ex , as one can see in Fig. \5\ c, d 

and f. The radiative decay time of the optically active excitons is inversely proportional to 

K ex and is directly proportional to the characteristic exciton decay time, as can be described 

by 

_1 _ 8e 2 um r P 2 _ V 

r ~ 9m 2 c 2 h To*' { } 

where n r is the refractive index of the media, c is the speed of light, P = —i(S\p z \Z) is the 
Kane transition matrix element, u is the light frequency, and T> is the depolarization factor, 
which takes into account the decrease of the magnitude of the light electric field when it 
penetrates into the NC. In CdSe structures Tq should be on the order of 1.7 ns.— This value 
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FIG. 5. Dependence of the square of the overlap integral, K , on the conduction band offset, Uos 
(a), and (e); and CdSe core radius, a (b), (d), and (f), calculated for an exciton (a) and (b); a 
positively charged trion and (d); and a negatively charged trion (e) and (f). 

was estimated using n r = 1.5, hu = 2 meV and 2P 2 /mo = 19.0 eV. The depolarization 
factor T> generally depends on the NC shape and structure. In giant CdSe/CdS NCs the 
depolarization factor can be estimated using V = [3e/(2e + ecds)] 2 for spherical CdS NCs, 
where e CdS = 5.35 is the high frequency dielectric constant of CdS and e = nl = 2.25 is 
the dielectric constant of the media surrounding the NC. This results in T> « 0.47 . 

The radiative decay times of the Bright band-edge excitons has a complicated dependence 
on the exciton exchange energy r\\ and the hole level splitting A .— The expression for the 
decay time is significantly simplified when A ^> rji . In this limit the radiative decay time, 
T ex , of the first optically active exciton can be written as 

1 K 

— = |f ( 49 ) 

'ex ZT 

This case is realized in giant CdSe/ CdS NCs, due to a small exchange energy. From Fig 15] one 
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can see that the overlap integral squared decreases by a factor of 2-3 when CBO decreases 
from 50 meV to 300 meV. This decrease of K ex is consistent with observed changes in the 
radiative decay time. One can see also from FigJS^ that an increase of the core radius leads 
to a strong localization of an electron in the core and increases K ex significantly. This in 
turn should shorten the radiative decay time. 

B. Positively charged trion transition probabilities 

The fine structure of the band-edge positive charge trions in NCs made of semiconductors 
with four-fold degenerate or quasi-degenerate valence bands (like CdSe, CdS, InAs, CdTe 
and etc.) generally consists of four doubly degenerate states Ey 2 , E 3 / 2 and Ei/ 2 and one 
four-fold degenerated state E$n , whose energies we described above. All of these states are 
optically active and can decay radiatively with different decay times. This in principle could 
lead to a very complex dependence of the decay time on the temperature as a result of the 
thermal population of the different levels within the trion fine structure. 

The resonant optical excitations of the positively charged trion in NCs with an extra 
hole and its radiative decay are determined by the probabilities of the optical transitions 
between one of a hole state, IS3/2 , with M = ±1/2, ±3/2 and the one of the trion states, 
Tr [±1/2^ (3/2, ±3/2); (5/2, ±3/2) ; (5/2, ±5/2) , and (1/2, ±1/2)]. These probabilities 
are proportional to the following sum of the matrix element squared: 

T(M,Tr) = \(^ h M(r hi )5(r e - r hl )\ep\^ g (r e ,r hl ,r h2 ))\ 2 

i,j = 1,2; i ^ j 

= 2|(^ / (r M )5(r e -r, 2 )|ep|$^(r e ,r w ,r ?i2 ))| 2 . (50) 

The sum takes into account the probability of an electron recombination with two holes 
( j = 1,2). The sum can be replaced by a factor of two because the matrix elements of the 
transitions with different hole states are identical, due the symmetrization of the two-hole 
wave function. Here e is the polarization of the emitted or absorbed light, and the momen- 
tum operator p acts only on the valence band Bloch functions. The calculation of T(M, Tr) 
with the wave functions is straightforward but cumbersome, because the hexagonal axes of a 
NC ensemble are randomly oriented relative to the light propagation of emission directions. 
First, let us consider transition probabilities T(M,Tr) 's responsible for the radiative decay 
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of the ground and first excited trion state. 

We will be considering here only the effect of linearly polarized light, which is sufficient 
for the description of the radiative decay. For linearly polarized light, ep in Eq. (150 j) can 
be written as 

ep = e z p z + - [e_p + + e+j5_] (51) 

where z is the direction of the hexagonal axis of the NC, e± = e x ± ie y , p± = p x ± ip y , and 
e Xj y and p XtV are the components of the polarization vector and the momentum operator, 
respectively, that are perpendicular to the NC hexagonal axis. 

In studied CdSe/CdS core/shell NCs, the crystal field and shape asymmetry parameter 
A is much larger than the trion exchange parameter rj . As a result the trion ground state is 
l/2~ . Let us calculate the transition probabilities T(3/2, l/2~; ir) and T(— 3/2, — 1/2~; ir) 
between the trion ground state and the hole ground state \M\ = 3/2 stimulated by the 
linear polarized light. Substituting Eq. (15 ip into Eq. (150]) and using the trion and hole 
wave-functions, after some straightforward but cumbersome calculations gives the following: 

T(-3/2, -1/2-; tt) = T(3/2, 1/2"; ir) = (S 2 ~ + B^f^^- sin 2 9 , (52) 

where 9 is the angle between the vector polarization of light e and the hexagonal axis 
of the NC, and where the trion overlap integral squared K tr + = \J drr 2 R^ =2 (r)i? ( r ) | 2 is 
determined now by the electron wave function i?f =2 of the trion, which is calculated for 
the Coulomb potential created by two holes. This radial function is described by Eq. (|29l) . 
with parameters a and £ that maximize the trion binding energy. The dependence of the 
overlap integral squared calculated for the positively charged trion is shown as functions of 
the CdSe core radius and CBO in Figs. 5 c and d. 

Equation (152"|) describes the probability of the trion ground state radiative decay into the 
ground hole state. There is nothing, however, that prevents the trion decay into the light 
hole excited states with M — ±1/2 . The photon emitted during such decay will be smaller 
(by an energy of about A ) than the photon emitted during the ground-to-ground transition. 
This energy will be absorbed by the hole. The probability of observing such a red "shake 
up" line is determined by the four independent transition probabilities: 

T(-l/2, -1/2-; tt) = T(l/2, 1/2"; tt) = (B 2 - B^) 2 ^fl cos * g , 

T(l/2, -1/2- tt) = IX-1/2, 1/2- ; tt) = (B 2 - B f^^- sin 2 9 , (53) 
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which are much smaller than the probability of the observation of the main line, because 
(B 2 - - B^f « (fl 2 - + B f for large A . 

Summing up the probabilities of the trion radiative decay into the ground and excited 
hole states over all light polarizations, we obtain the radiative recombination rate of the 
trion ground state as 

1 Ktr+ (54) 



2r ' 

One can show that the uppermost trion excited state l/2 + has the same radiative decay 
rate as the ground one ( Ty 2+ = t^I- ) and they are different from the decay time of the 
first optically active exciton state only by the overlap integrals (see Eq. fj4*9|) ). Comparison 
of the overlap integral squared for excitons and positively charged trions in Fig. E] shows 
that in NCs with small conduction band offsets they became much larger for the trions due 
to stronger attraction of electrons to the CdSe core. Henceforth in this paper we will use 
notations for the decay times, which are related to the notation of corresponding levels of 
the positively charged trion: t^1± , an d rfj^ an d r£j 2 . 

We find the probability of the transitions from the excited trion states (5/2, ±5/2) and 
(5/2, ±3/2) to the hole M = ±3/2 state to be 

T(±3/2, (5/2, ±5/2); vr) = Ktr+P% sin 2 9 , T(±3/2, (5/2, ±3/2); tt) = 2Ktr+p2 cos 2 g . (55) 

6 3 

The transition to the hole ground state has an energy of about A larger then the main 
transition from the ground state. 

The total transition probability from the trion ±5/2 states to the hole ground ±1/2 
state is described by the sum of the two probabilities: 

T(l/2, (5/2, 3/2); tt) = T(-l/2, (5/2, -3/2); vr)) = sin 2 9 . (56) 

In the same limit of large A , transitions from the trion excited state to the hole excited 
state have nearly the same energy as the main transition from the ground to ground states. 

When we average all possible crystal orientation and sum up Eqs. (1551) and (1561) . we can 
find the radiative recombination rate of the trion excited state E5/2 '■ 

1 3K tr+ (5?) 



7-5/2 5r ° 



Similar calculations give us the radiative rates for E1/2 and E 3 / 2 excited trion states: 

1 2K tr + 1 llK tr + 
3r r 3 ^2 



(58) 

r 1/9 6t T, /9 LbTQ 
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The situation is significantly simplified in studied here CdSe/CdS core/shell NCs. Due to 
small value of fj < A, the second trion level is eightfold degenerate with an average decay 
rate of 

Ktr+ (59) 




r 2 tr+ 4 \ r *C ' tK ' tX ) 2r 



which coincides with the decay rate of the uppermost and lowest trion states 1 / t*^ 2± (see 
Eq. ( 154]) ). As a result, the positive trion radiative decay time in CdSe/CdS core/shell 
NCs does not depend on thermo-population of the fine level structure and is temperature 
independent. 



C. Negatively charged trion transition probabilities 

The radiative decay time of the negatively charged exciton leaves a resident electron 
in NCs. In CdSe/CdS core/shell NCs, the wave function of the resident electron, R? e , is 
very different from the electron wave functions of the negatively charged trion, which is 
mainly controlled by the Coulomb interactions. The different size of the electron and hole 
wave functions significantly decreases the overlap integral, which controls the rate of trion 
radiative decay. 

The radiative decay and the resonant optical excitation of the negatively charged trion 
are determined by the probabilities of the optical transitions between the resident "free" 
electron state, with spin projection S z = ±1/2 and radial wave function R® , and one of 
the trion states characterized by the hole momentum projection M = ±3/2, ±1/2 . These 
probabilities are proportional to 

T(S Z ,M) = 2|(^(r 1 )5(r 1 -r e2 )|ep|^(^)* 2e (r el ,r e2 ))| 2 , (60) 

where the factor 2 is due to the sum over electron coordinates. The ground negative trion 
state is formed by the hole with M = ±3/2 . The probability from the ground state is 

T(-l/2, -3/2; tt) = T(l/2, 3/2; tt) = sin 2 9 . (61) 

Here, 

Ktr- = \ [ [ r 2 1 dr 1 r 2 2 dr 2 R° e (r 1 )R 2e (r 1 ,r 2 )R Q (r 2 )\ 2 , (62) 
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where i?2e is the symmetrized two-electron wave function in Eg. (I44p . Rq is the wave func- 
tions of a single hole in Eq.fl2]), and R® is the wave function of a single resident electron in 
Eq.([I7D. 

Figures 5e and 5f show the dependence of the overlap integral squared, K tr - , for the nega- 
tively charged trion on the CdSe core radius, a , and the CBO, | C/off I ■ One can see from these 
figures that K tr - always increases with a and \Uo&\ ■ At the same time the figures show 
that for each a there is critical |£/off| an d vise versa where K tr - vanishes. These critical 
points are connected with the derealization of the resident electron into the CdS shell. As we 



show above, the derealization condition is described as Wq v = a^2mcdSc\Uos / h ~ 1.296. 
At this condition the spread of the electron wave function into the CdS shell vanishes the 
overlap integral. 

Taking a sum over all light polarizations in Eq. flST]) . we obtain the radiative decay rate 
from the negative trion ground state as 



1 _Ktr- 



(63) 



As in the case of the positively charged trion, the radiative decay time of the negatively 
charged trion is different from the radiative decay time of the first optically active exciton 
only by the overlap integral. For the negative trion excited states with M = ±1/2 , we have 
following transition probabilities: 

2K -P 2 

T(-l/2, -1/2; tt) = T(l/2, 1/2; vr) = cos 2 9 , (64) 

T(l/2, -1/2; vr) = T(-l/2, 1/2; tt) = sin 2 9 , (65) 

6 

Taking a sum of the transition probabilities over all states and light polarizations, we obtain 
the radiative decay rate of the negative trion with M = ±1/2 : 

T l/2 2r 

The radiative decay time of the trion excited state, Ty 2 , is the same as one for the ground 
state: Ty 2 = r^ 2 . 

IX. DISCUSSION 

This paper has developed a theoretical background for qualitative and quantitative anal- 
yses of multiple experimental measurements conducted in CdS/CdS core/thick shell NCs 
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and CdSe/CdS dots-in-rods nanostructures (NS). The major unusual properties of these NS 
are connected with the relatively small, and temperature and strain dependent conduction 
band offset, which provides only weak confinement of electrons in the CdSe core, while the 
large valence band offset in these NS leads to the strong confinement of holes. As a result, 
the Coulomb potential created by strongly confined holes plays an important role in electron 
confinement. 

The optical properties of such NS should be considered to be in the intermediate con- 
finement regime.— The calculation procedure consists of several steps: (i) finding the en- 
ergy spectra and wave function of strongly confined holes; (ii) calculation of the adiabatic 
Coulomb potential creating by the hole charge distribution; and, finally, (iii) calculation the 
energy spectra of the electron moving into the total potential created by the conduction 
band offset and by the Coulomb potential. In this paper, we realized this program assuming 
that the thickness of CdS shell is much larger than the radius of electron localization around 
the CdSe core. The limit is satisfied in "giant" CdSe/CdS core/shell NCs grown by the 
Los Alamos^ - - 1 *^ and Deburter-^ 1 ^ groups. The theory could be extended, however, to 
quantitative consideration of CdSe/CdS dots-in-rods NS and/or to take into account the 
finite thickness of the CdS shell. 

Some effects of the finite CdS shell can be taken into account quantitatively on the basis 
of the current calculations if the radius of electron localization around the CdSe core, (a , 
is smaller than the CdS layer thickness. For example, the negative trion ionization energy 
should be increased on the energy of the first confined level in CdS NCs with the same shape 
and size as original CdSe/CdS core/shell NCs. The radiative decay of the negative charged 
trion also does not vanish in NC with finite shell thickness when a derealization condition 
is reached. The time can be calculated by replacing R® in Eq. fl62l) with the wave function 
of the the first confined level in CdS NCs. 

In the paper we have also derived general expressions for the radiative decay time of the 
band-edge trions in various NCs. The calculations show that these times are proportional to 
2ro/K tr ± . The major difference in these times is connected with the square of the overlap 
integral between electron and hole wave functions in positively charged K tr + and negatively 
charged K tr - trions. In CdSe/CdS core/thick shell NCs where the electron exchange energy 
becomes smaller than the temperature, the radiative decay time of excitons is two times 
longer than 2r / K ex given by Eq. (fl9"j) due to the equal population of the Dark and Bright 
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exciton states. 

A detailed comparison of our calculations with experimental data on the optical properties 
of the CdSe/CdS core/shell NCs requires a knowledge of the temperature dependence of 



14 



is not 



the conduction band offset. This temperature dependence suggested first in Ref. 
universal, because it has a strain contribution^ and should depend on the sample preparation 
technique. Nevertheless, the decrease of the radiative decay time by a factor of two observed 



in Refs. B and 



9 
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is consistent with an increase of the overlap integral squared with the 
CBO increased from 100 to 300 meV. The small splitting of the Dark/Bright exciton states 
measured in the giant CdSe/CdS NCs^ is also in qualitative agreement with our calculations. 

In summary, we have developed a theory of the optical properties of CdSe/CdS core/thick 
shell NCs. The calculations takes into account the complex structure of the CdSe valence 
band and inter-particle Coulomb and exchange interaction. The theory describes the fine 
structures and radiative decay time of excitons and positively and negatively charged trions. 
The results of our calculations are in qualitative agreement with experimental data available 
for these NCs. We believe that calculated energy spectra of negatively and positively-charged 
trions will explain suppression of nonradiative Auger processes observed in these structures. 
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Appendix: Normalization constant 

From the normalization condition for an electron bound to the singly or doubly charged 
core we obtain the following analytical expression for the normalization constant squared: 

B 2 = 256 [2 (a 2 -l)C + aV] 2 

a 3 [q (/i, a, £) + v^vrexp (2/a 2 ) erf [\f2/a) p (fx, a, Q\ 
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where p(fi,a, () and q (//, a, () are given by 



q ( Atj a , C) = 16«C 2 (4 - 55a 2 ) + 256c< 3 (2 + 2( + ( 2 ) + 16a 3 (fi (27a 2 - 4) 

+ 4a V (4 - 15a 2 ) 
p a , C) = 4C 2 (16 - 56a 2 + 55a 4 ) + 4a 2 C/i (40a 2 - 27a 4 - 16) 

+ aV (16 - 24a 2 + 15a 4 ) 



and erf (x) is the error function defined as 



2 



X 



erf (x) 



e~ t2 dt 
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